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Abstract 

The influence of disordering upon critical behavior of the system with hidden degrees 
of freedom is considered. It is shown that there is a tricritical behavior in the 
constrained system, while in the unconstrained system only phase transitions of the 
second order occur. 
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1 Introduction 



Recently, phase transitions under a constraint imposed on some "hidden" de- 
grees of freedom coupled to the critical ones have again attracted attention 
[1]. As it is well known, many real systems possess a hidden degree of freedom 
[2-4], which can be described as the nonfluctuating parameter in the vicinity 
of the phase transition. Fisher [5] stressed the importance of the constraint 
upon the hidden degree of freedom. Fisher's theory of phase transition in con- 
strained systems was generalized [6,7] to include the possibility of first- order 
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phase transitions. In particular, it was found that the tricritical behavior co- 
incides with the critical behavior of the ideal system provided a > (a is 
the critical exponent of the specific heat) and with the renormalized one pro- 
vided a < 0. Early renormalization group studies of systems with constrained 
nonfluctuating parameter were presented in Refs. [4,8,9]. 

As it is shown in Ref. [1], the nature of a phase transition may strongly depend 
on the constraint imposed on the system. For the case of the unconstrained 
system the stability criterion for the occurrence of continuous phase transition 
can be found within the mean field theory. On the other hand, for the case of a 
constrained system a different criterion, which is distinct from the mean-field 
one, is calculated. This criterion results in a fluctuation-induced renormalized 
Heisenberg tricritical point. 

The main purpose of this paper is to consider the influence of disordering upon 
critical behavior of a system with hidden degrees of freedom and to focus at- 
tention on a criterion for the occurrence of continuous phase transition in this 
system. Such disordering can be caused, e.g., by the presence of "frosen" impu- 
rities [10]. The analysis of renormalization group equations for random system 
with constrained nonfluctuating parameter has been presented in Ref. [11]. It 
was found the fixed points and shortly described the critical behavior. Here 
we discuss the tricritical behavior and especially a first-order transition taking 
into account results of [1]. 



2 Renormalization group equations 



Let us consider the disordered system in which an n-component vector or- 
der parameter S(x) is coupled with the scalar nonfluctuating order parameter 
2/(x). This system can be described in momentum space by the effective Hamil- 
tonian [11] 



^ = 2/ (2^(2^ r(qi ' q2)S(qi)S(q2) + 

+ / (S? M(qi ' q2 ' q3 ' ^) S Cq 1 )S(q 2 )S(c 1 3)S(q 4 ) + 

+ / W^WfWf^ q2) c i3)?/(qi)s(q 2 )s(q3) + 

4 / ^yWf^ q 2 )y(qi)y(q 2 ) + / ^Mq)i/(q)- (i) 

The Hamiltonian (1) can be used to study the ferromagnetic phase transition 
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in the three-dimensional Hubbard model with frosen impurities. Herewith, the 
coupling p is pure imaginary 

li = i(kT e ) 1 ' 2 U*' 2 N'(E f ), (2) 

where N(Ef) is the density of states on the Fermi energy Ef and U is the 
Coulomb potential. The study of this case is of particular interest. 

Following the standard method of the renormalization group [12] and its exten- 
sion to random systems [10] the recurrence relations for the potential averages 
over the probability distribution function can be constructed by averaging 
on the recurrence relations for the potentials of the Hamiltonian (1) of the 
inhomogenious system and for the second cumulants: 



r' = b 2 {r + [4(n + 2)u - A + 2nz}A(r)}, (3) 

u' = ¥{u - K 4 In 6[4(n + 8)u 2 - QuA}}, (4) 

A' = b e {A - K 4 In b[8(n + 2)uA - 4A 2 ]}, (5) 

z =¥{z - K 4 \nb[8(n + 2)uz + 2nz 2 -2zA]}, (6) 

w' = ¥{w - K 4 In 6[8(n + 2)uw + 2nw 2 + Azw - 2wA]}, (7) 



where we introduce the following symbols for second cumulants: A ~< rr > c , 
v ~< rh > c , p ~< hh > c , and use the notation: u = u—^, A = A—^ + ^j^, 

2 2 2 

z = ^- — w = All other definitions are standard in renormalization 

P PO ' PO 

group theory: = ^ is a quality proportional to the surface area of the 
unit sphere in d = 4 space, b is the scale parameter and A(r) is the integral 
of a closed loop of pair correlation functions (see, also, [13]). 

In Eqs.(3-7) we separate the coefficient of the nonfluctuating parameter y(q = 
0) from those of y(q ^ 0) because of their possible role in constraining systems 
[4,8]. We make a shift of the variable y(q_ = 0) for the disappearance of the 
linear term in the variable y on each step of the renormalization procedure. 
It is also necessary to write the equations for the cumulants v, u , p/i, poPo- 
However, the structure of the new equations is such that they do not violate the 
stability of fixed points, which were defined by Eqs.(3-7). Using the definition 
of cumulants A, v, p one can show that A > 0. 

Before discussing the renormalization group analysis, consider the mean-field 
theory results. After integrating over the nonfluctuating parameter in (1), we 
obtain a new effective Hamiltonian for order parameter S(q). A solution of this 
new integrated out effective Hamiltonian in the mean-field approximation gives 
the boundary conditions of instability for the Hubbard model: \ s — \^ = 0, 
A^°) = 0, where \ s = u and Aji ) = //q/2/3o, for unconstrained system (A^ = 
A^ = A) and A s = 0, A^ = 0, where A^ = p 2 /2(3, for the constrained system 
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Fig. 1. Phase diagram in the mean-field approximation for the unconstrained system 
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Fig. 2. Phase diagram in the mean-field approximation for the constrained system 

(A(°) = 0). These conditions divide the first and second order phase transition 
ranges in (A c , A^)-A s -plane (Fig. 1 and Fig. 2, see, also, [1]). Using the Stoner 
criterion r s = 1 — (U/2)N(Ef) = 0, we have the condition for a continuous 
phase transition for the unconstrained system 



while the first order phase transition occurs for the contrary case. For the 
constrained system we have the continuous phase transition, if 



According to the mean-field theory for random systems, the impurities do not 
influence on the stability ranges of the system. 




(8) 



N"(E f ) < 0. 



(9) 



A feature of the recurrence relations is a closed system of two equations: (4) 
and (5). It is easy to find the fixed points of this system. We have as usual 
the Gaussian fixed point (G) with u* = A* = 0, the Heisenberg fixed point 




Fig. 4. Flow trajectories in the A = A* plane 

(H) with u* = e/AK^in + 8), A* = 0, the non-physical fixed point with 
u* = 0, A* = — e/4K 4 , and random fixed point (R) with u* = e/lQK A (n — 1), 
A* = e(4 — n)/ 8X4(11— 1). Using these fixed points, one can easy find the rest 
fixed points from the recurrence relations. 

The full list of fixed points, the critical behavior, which is determined by their 
values, and their stability have been presented in [11]. We can see from the 
Table in [11] that some parameters and their selfvalues Aj change sign at n — 4. 
This changing in sign strongly influences the forms of flow trajectories, and 
we are obliged further to discuss separately the n < 4 and n > 4 cases. It is 
convenient to use the two parameter spaces: (A s , A c , A) for the unconstrained 
case (A^ -* = A^, 1 -* = A c ) and (A S ,A^,A) for the constrained case (A^ = 0). 
Moreover, it is turn out that the interesting fixed points are in the planes A = 
and A = A* 7^ in both spaces. Due to the relation u = u — fi 2 /2(3 = A s — A^ 
all fixed points align on two parallel lines in these planes. The flow trajectories 
are identical in both spaces, so we can demonstrate them simultaneously in the 
A = (Fig. 3) and A = A* (Fig. 4) planes in the spaces of the renormalized 
coupling constants g s , gf\ and The notation of fixed points, which 

correspond to the subspace (A s , A^, A) are written in brackets. It is important 
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to point out that fixed points with z* = 2 (\^* — A^*) determine the 
renormalized critical behavior in accordance with [5]. 

One can see from Fig. 3 and Fig. 4 that the most stable fixed points and, 
hence, the critical behavior is independent of the existence of constraint. One 
point is H 2 , when n > 4, and the critical behavior is ideal for a pure system 
and another point is R 2 , when n < 4, and the critical behavior is determined 
by the presence of the frosen inhomogeneities. Thus, the behavior coincides 
with the one found in [10]. 

The critical behavior at the boundary of the stability region is determined by 
the fixed points R\ (for n < 4) and Hi (for n > 4) in the system without 
constraints, and by the "renormalized" fixed points RR\ (for n < 4) and RH\ 
(for n > 4) in system with constraints. It is easy to see that the boundary of 
the stability region in constrained systems reduces the criterion (9) for pure 
system to 

3(^(£/)) 2 OLH/nVH „ 

2N(E f ) < e/(n + S) + a H /nu H N (Ef) + °^ 



where an = 2 (n+8) 6 an< ^ v h — \ J r 4(^+8) 6 are ^ ne specific heat and the cor- 
relation length critical exponents for the Heisenberg fixed point, respectively 
(see, also, [1]). By analogy with a pure system, one can find the corresponding 
criterion for a random constrained system in the form 

2N(E f ) < e/A(n-l) + a R /nu R N ™ + °^ (U) 



where a R = 8 ™~^ e and = \ + 32 (n-i) e are corresponding exponents for 
the random fixed point. Thus, in the constrained system the flow trajectories 
run away in both the pure system fixed point Hi and the random fixed point 
Ri, and, hence, the critical behavior at the boundary of the stability region 
reduces to the tricritical behavior. 



3 The tricritical behavior 



Using the standard field theory and renormalization group method (see, for 
example, [14]), we obtained the next equation for the free energy in the one 
loop approximation 



(3F(T, M) = \r s {T)M 2 + i(A, - X^)M' - 
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2 J r s (T) + q 2 4! 2 J (r s (T) + g 2 ) 2 + 
1/ \ f i if XiM 2 /6 x 

1 /■ , , /„ A 2 M 2 /6 x 

where M = \ J2 a M 2 ) is the magnetization, and 

Ai = i(A.-A?>), A 2 = i(3A s -Af-2A«). (13) 



Note that here we used coefficients of Eq.(l) in form: u = A s /4!, ^q/2(3q = 
A^/4!, /i 2 /2/5 = A^/4!. To get the renormalized free energy, we used the 
theory of Refs.[14] and [1]. Finally, we have 



(3F(t, M) - (3F(t, 0) = -tM 2 + ^(~9s- ^ 0) )M 4 + 

(n - l)f(t + \ gi M 2 ) + fit + \g 2 M 2 ), (14) 

where t, g s , and are the renormalized by frozen impurities values, g\ and 
g 2 are renormalized coupling constants Ai and A 2 and 

f{x) = \x 2 {\nx- 1 -). (15) 

For the unconstrained case g^ = g^ = g c we have the conventional 4 -model 
with an effective coupling constant proportional to g s — g c and the runaway 
flow trajectories do not intersect the boundary of the stability range. Hence, we 
can suggest that the phase transition remains second order. It should be noted 
here, that the fixed points Hi and R\ are not tricritical points. In particular, 
this fact is in accordance with the result for pure systems [1]. Also, following 
[1], we can look on the lines connecting the origin and fixed points Hi and Ri 
in Figs. 3 and 4 as boundary lines between two second order phase transition 
regions with different critical behaviors. 

For the constrained case gf^ ^ g^ the runaway flow trajectories intersect 
the boundary of stability range. The critical behavior at the boundary of that 
stability region is determined by the points RRi (n < 4) and RHi (n > 4). We 
suppose that this means rather the smeared phase transition than the phase 
transition of first order, in accordance with [10] and [15]. 
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In summary, we considered the influence of disordering upon the critical be- 
havior in systems with hidden degrees of freedom taking into account the 
possible role of the constraint upon these degrees of freedom. It is shown that 
there are the tricritical fixed points in the constrained systems and we have 
rather smeared phase transition than the first order transition, while in the un- 
constrained system all phase transitions are the phase transition of the second 
order. 

References 

[I] M.L.Lyra, M.D.Coutinho-Filho, and A.M.Nemirovsky, Phys. Rev. B 48 (1993) 
3755-3769. 

[2] D.R.Nelson and M.E.Fisher, Phys. Rev. B 11 (1975) 1030. 
[3] J.M.Kosterlitz, D.R.Nelson, and M.E.Fisher, Phys. Rev. B 13 (1976) 412. 
[4] C.M.Chaves, P.Lederer, and A. A. Gomes, J. Phys. C 10 (1977) 3367. 
[5] M.E.Fisher, Phys. Rev. 176 (1968) 257. 

[6] Y.Imry, O.Entin-Wohlman, and D.J.Berdman, J. Phys. C 6 (1973) 2846. 

[7] V.Dohm, J. Phys. C 7 (1974) L174. 

[8] Y.Achiam and Y.Imry, Phys. Rev. B 12 (1975) 2768. 

[9] Y.Achiam, J. Phys. C 10 (1977) 1491. 

[10] T.C.Lubensky, Phys. Rev. B 11 (1975) 3573. 

[II] V.M.Laptev and Yu.N.Skryabin, phys. stat. sol. (b) 91 (1979) K143. 
[12] K.G.Wilson and J.Kogut, Physics Rep. C 12 (1979) 75. 

[13] Yu.A.Izyumov and Yu.N.Skryabin, Statistical Mechanics of Magnetically 
Ordered Systems (Consultants Bureau, New York, London, 1988). 

[14] D.J.Amit, Field Theory, The Renormalization Group and Critical Phenomena 
(World Scientific, Singapore, 1984). 

[15] A.Aharony, Phys. Rev. B 12 (1975) 1038. 



8 



